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Abstract. Recent lattice-QCD results for the dressed-gluon propagator are used within the quark 
£Nl , Dyson-Schwinger equation to determine the gluon-quark vertex dressing necessary to reproduce 

' the lattice-QCD results for the dressed-quark propagator. Both quenched and full QCD lattice 

simulations, for a range of low quark current masses, are analyzed. The chiral extrapolation is made 
y— 5 . through this continuum DSE form. Resulting chiral and physical pion observables are investigated. 
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c3 : INTRODUCTION 

O ; 

Lattice-QCD simulations address hadronic observables via matrix elements of currents 
and sources; necessary approximations and truncations introduce systematic errors from: 
a finite volume of discretized space-time, unphysically large current quark masses, and 
43 the quenched approximation. While steady progress is being made in the reduction 

of such errors, it is useful to gain insight into dominant field theory mechanisms by 
comparing to hadronic results from covariant modeling of continuum QCD. The Dyson- 
Schwinger equations (DSEs), the equations of motion of the theory, provide such an 
opportunity. Here systematic errors can arise from truncations that replace high-order 
correlations by infrared phenomenology of low-order n-point functions. Symmetries can 
^ ■ provide significant quality control. The dressed quark and gluon propagators and low- 

order vertices of the theory are important elements of the kernels needed for the bound 
state equations: Bethe-Salpeter equation (BSE) for mesons and Faddeev equation for 
baryons. Prior to the last few years, one had only the known ultra-violet behavior of the 
n-point functions as a guide. 

Recently, full-QCD lattice results for gluon [1] and quark [2] propagators have been 
made available; these are two of the three nonperturbative quantities linked by the 
quark propagator DSE. Phenomenological information on the third quantity, the dressed 
quark-gluon vertex T^(k,p), can therefore be inferred. We have previously carried 
out such an analysis of quenched-QCD lattice propagators [3] and the quark-gluon 
vertex [4]; the procedures and technical details employed here are described in Ref. [3]. 
Here, and in the lattice simulations, Landau gauge and the Euclidean metric is used. 
Briefly, the approach is the following. From the renormalized quark DSE, we have 
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FIGURE 1. Left Panel: Full-QCD lattice results [2] reproduced from the full-QCD lattice gluon prop- 
agator [1] via the continuum DSE. Note vertical log scale. Values of lattice m and ma are shown; Right 
Panel: Fourier transform of the Dirac scalar part of the chiral S(p). Cusps indicate confinement. 

=^2iy- P+ Z4m(fl) +£'(/?, A) and the regulated self-energy is 

Z'(p,A)=Z 1 J^g 2 D^(p-q) ^ Yll S{q)T a v {p-q,p) . (1) 

Here A is the regularization mass scale, the Z;(/i 2 ,A 2 ) are the usual constants 
from renormalization at scale /l. Two are required so that as p 2 — > jj, 2 , we 
have S~ 1 (p) — > iy ■ p -\- m(ii) . The lattice results for S(p) are given in terms of 
S(p) — Z(p 2 , ii 2 ) I (iy- p + M(p 2 )). We compare by solving (1) with the kernel factors 
Zi g 2 D^ v (k) T a v (Jfc, p) replaced by D l {* (k) 4£ y v H (k 2 ) . Here D l {* (k) is a fit to the lattice 
gluon propagator and the quantity Ti(k 2 ) is a phenomenological vertex amplitude de- 
termined so that the DSE solution for S(p) fits the lattice data. In general the transverse 
vertex has eight amplitudes and a dependence upon two momenta; however the data 
being fitted here do not warrant more. We note that one may identify the kernel as 
4na e ff(k 2 ) v (k) where is the th order gluon propagator. We ensure that the 
leading log behavior of all quantities conform to the 1-loop renormalization group 
behavior of QCD. 

RESULTS 



TABLE 1. Chiral condensate and pion decay constant (chiral 
f%, physical f K ) from the lattice-guided DSE kernel. 
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Full-QCD 


Qu-QCD 


4?) fi=l GeV 


-(0.24 GeV) 3 


-(0.23 GeV) 3 


-(0.19 GeV) 3 


fn 


0.090 GeV 


0.072 GeV 


0.063 GeV 


fx 


0.092 GeV 


0.075 GeV 


0.066 GeV 



Our fit to the full-QCD Df v {k) uses the "model A" form previously employed for 
quenched data [5] but now with Nf = 3 and the new parameter values: A = 3.25, 
A g = 0.54, a = 1.15 and Z g = 1.22. Our fit gave priority to M(p) for the four available 
values of lattice ma shown, along with the results, in the left panel of Fig. 1. The 
parameterized form used for Ti(k 2 ) is the same as we have previously used [3] for 
the quenched lattice case. The (dimensionless) parameters of Ti(k 2 ) found here for 
the full-QCD case are: a\ = 4.5, a2 = 2. 1,03 = 18.1,fr = 0.31. The result of the chiral 
extrapolation provided by the DSE kernel is also shown. The lattice current masses are 
equally spaced and, in the region p > 3GeV, the lattice results are only approximately 
so; the ma = 0.01 case deviates most from the pattern. The DSE fit has been made with 
the constraint that the M(p, m) approach the correct ratio. 

If a propagator of a field theory in Euclidean metric violates the Osterwalder-Schrader 
axiom of relection positivity [6], then this is a sufficient condition for confinement of the 
corresponding excitation [7]. In the right panel of Fig. 1 we display the magnitude of 
the Fourier transform of a s (p4,p = 0), the Dirac scalar amplitude of the chiral limit 
quark propagator. For a free particle with mass m, A$(T) <*exp(— mT). The cusps 
indicate changes of sign and thus confinement in both quenched and full-QCD. The 
dotted line corresponds to a propagator that has a single pair of complex conjugate poles 
at p 2 = -0.309 2 ±/0.192 2 GeV 2 . The dashed curve corresponds to the ladder-rainbow 
model that describes a large variety of light quark observables [8]. 

The form of the deduced DSE kernel allows a chiral symmetry-preserving Bethe- 
Salpeter kernel to be obtained as in the ladder-rainbow case [9]. This produces the 
chiral physics observables shown in Table 1 . Full-QCD evidently produces an acceptable 
condensate, while the quenched approximation underestimates by a factor of two [3]. 
The values of f n , both chiral and physical, are marginally improved by full-QCD but 
they remain about 15% too low. 
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